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ABSTRACT 

A model is presented which assumes the existence of cold dense clouds near the central engine of active 
galactic nuclei (AGNs). The effects of such clouds on the observed spectrum are explored. It is shown that this 
model is consistent with the complicated observed spectra and variability behavior of most extensively studied 
Seyfert nuclei. The results are compared with other proposed models. The existing observational evidence 
appears to support the “ cloud-model.” 

Subject headings: dense matter — galaxies: active — galaxies: nuclei — galaxies: Seyfert — quasars: general 


1. INTRODUCTION 

The generally power-law spectrum of a class of AGNs rep- 
resented by radio-quiet quasars and Seyfert nuclei is consistent 
with a synchro-Compton model where the primary synchro- 
tron radiation is modified through Compton scattering by a 
nonthermal electron-positron pair plasma near a supermassive 
black hole (e.g., Fabian et al. 1986). Recently additional fine 
structures, such as the Fe line-edge feature and hard X-ray 
bump, have been observed in many of Seyfert nuclei (e.g., 
Nandra et al. 1991 and references therein). These observed 
characteristics are naturally explained by the “ two- 
component” model which consists of the hot nonthermal 
plasma responsible for the Comptonization of the primary syn- 
chrotron radiation and the cool plasma which reprocesses and 
reflects the nonthermal radiation (e.g., Guilbert & Rees 1988, 
hereafter GR88; Lightman & White 1988, hereafter LW88). 
The possibility for the cold matter to lie near the central region 
is supported by the simultaneous rapid variations of the X-ray 
continuum and Fe line fluxes, to within ~300 s, recently 
detected in Seyfert nucleus NGC 6814 (Kunieda et al. 1990). 
The cold component is envisioned to be either in a form of a 
thin disk (or slab) or an assembly of clouds distributed over a 
three-dimensional configuration (LW88, GR88). For conve- 
nience, we shall call the former kind “the disk model” and the 
latter “ the cloud model.” In order for the cold clouds to exist 
near the central engine they must be sufficiently dense, and 
some physical pressure is required to confine them. This pres- 
sure may be supplied by the equipartition magnetic fields (Rees 
1987; GR88; Celotti, Fabian, & Rees 1992, hereafter CFR92). 

The spectral and temporal behavior of the outcoming radi- 
ation depends on the geometry of the radiating and refecting 
plasmas. So far, numerical calculations of reflection of non- 
thermal radiation from cold matter have been carried out only 
for the disk (or slab) geometry (e.g., LW88; George & Fabian 
1991). However, existing observations already appear to favor 
the cloud model (see § 4, also Tsuruta 1992). In this paper, 
therefore, we extend the previous work by carrying out 
numerical calculations of reflection of nonthermal synchro- 
Compton radiation from cool, dense clouds confined in a geo- 
metrically defined region. 

Section 2 outlines the basic equations and methods adopted 
in this work. The results are presented in § 3. Discussion and 
concluding remarks are given in § 4, the last section. 


2. BASIC EQUATIONS 

In deriving our equations the general approaches of GR88, 
LW88, and CFR92 are followed. Further details are found in 
Sivron & Tsuruta 1992 (hereafter ST92). The free parameters 
used in our work are the total bolometric luminosity L, the 
mass of the black hole M, and the total covering factor due to 
the cloud system/ 

In this work we consider accretion-powered AGNs radiating 
close to the Eddington limit, which are generally thought to 
apply to Seyfert nuclei and radio-quiet quasars. The total 
luminosity is therefore L = qtitC 2 , where r\ is the efficiency 
of the process and tit is the accretion rate. For a continuously 
smooth and isotropic i nflow of matter the accretion rate is 
tit = 4nR 2 nm H (v ac Jv {{ )92GM/R. The density of the infalling 
matter is ~m H n(R), where n is the average number density and 
m H is the mass of one hydrogen atom. (Assuming cosmic abun- 
dance this is a reasonable approximation for the mass of an 
atom.) The infall velocity of the accreted matter is u acc , which is 
assumed to be proportional to v„, the free-fall velocity. Hence 
the minimal optical depth of isotropically infalling matter is 
determined by r = J k cs nm H dR , which can be expressed as 



where R g — IGMjc 1 is the gravitational radius of the central 
black hole, and t acc and t ({ are the accretion and free-fall time 
scales, respectively (GR88). The Eddington luminosity is 
•f'Edd — 47i GMm H c/a T , where G, c, and c T are the gravitational 
constant, the speed of light, and the Thomson cross section. 
When the source is radiating close to L Edd , therefore, x will 
exceed one near the center, and the original power-law emis- 
sion will be shrouded and reprocessed by the infalling matter. 
Consequently the final spectrum is expected to deviate from a 
simple power law (GR88, LW88), as indeed observed by the 
X-ray satellites such as Ginga (Nandra et al. 1991, and refer- 
ences therein). 

When the inflow of matter is not smooth the above equa- 
tions should be modified. In the two-component cloud model 
the frequency dependent optical depth of one cloud is 

[•K + 2J c ,d 

T vicid = I dR(K v + K KN )n lh m H , (lb) 
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where n th m H is the density of the cloud, k kn and k v are, respec- 
tively, the free-free and Klein-Nishina opacity, v is the fre- 
quency, and / cld is the radius of a cloud. An accurate optical 
depth that includes the atomic properties and the Einstein 
coefficients, such as the one presented by Ferland & Rees 
(1988), will be included in our subsequent work (ST92). 

Because of fast Compton and faster bremsstrahlung cooling 
most of the shrouding infalling clouds should be relatively 
cold, but the temperature of the clouds cannot be lower than 


Z M = 


(1-^4) 


dv 


16a„ nR 2 


11/4 

(1 _ e -Tvu.d ) e ~M) (2) 


because of the absorption of the primary radiation and the 
subsequent thermal reemission. Here A is the albedo (per 
cloud), cr B is the Stefan-Boltzmann constant, and f(R) is the 
covering factor at R, which might be also frequency-dependent. 
This temperature is reached if n, h , the number density of matter 
in the cold clouds, is high enough so that the bremsstrahlung 
cooling time satisfies t hTemss <$ t Compton <$ t hIolton , where 
f Compton an d Synchrotron are the cooling times due to Compton 
scattering and synchrotron emission, respectively. This high 
density can result from a confining pressure supplied by the 
equipartition magnetic fields, P magnetic = B 2 /16n ~ nm n GM/R 
(Rees 1987). Then a two-phase structure with dense cold clouds 
and hot nonthermal plasma can exist. By equating the gas 
pressure to the confining magnetic pressure we get (GR88, also 
see LW88) 


if =mHc2 G8 ■ e> 

Another important relation between the above variables can 
be derived from arguments similar to those that lead to equa- 
tion (1). The accretion rate of the infalling clouds is the product 
°f n c\d> the number of clouds per unit volume, the mass of one 
cloud, the accretion velocity of the clouds, and the surface area 
of a spherical shell. When the clouds hold most of the infalling 
mass, their accretion velocity is ~ v acc so that 


M 


dM cld 

dt 



4nR 2 n ( . id n th m H v acc . 


( 4 ) 


One of the assumptions made for the cloud model is that a 
significant fraction of M does not fall into the equatorial plane 
to form an accretion disk when the clouds are small and dense 
enough, even when a small amount of angular momentum is 
present. Therefore, when the time scale for the cold matter to 
fall into a disk is longer than the typical infall time scale R/v ff , 
most of the matter will be in the form of clouds. We assume 
that the typical collision time f coll < l/(n cld nl 2 Id v cld ) is equal to 
the time scale for falling into a disk, where v cJd is the random 
velocity of the clouds. From this requirement we obtain the 
following relation between n cld and / cld . 


‘cld ‘cld 



If the main force acting on the clouds is gravitational, u c!d is 
roughly the Keplerian velocity. Otherwise the explicit 
dynamics should be considered. 

As the input to our spectral calculations, a point source is 
adopted which emits a power-law spectrum from the inner- 
most region of the central engine. The model has to satisfy the 
variability requirement imposed by the absence of IR-optical 


variations during violent X-ray flares, which is noted for 
Seyfert nucleus NGC 4051 (Done et al. 1990). Therefore, we 
assume that the IR-optical portion of this primary radiation is 
effectively absorbed by the extremely small cloudlets in the 
vicinity of the black hole which are transparent to the higher 
energy photons, as suggested by CFR92. We shall call these 
small clouds the CFR cloudlets.” In our model this modified 
spectrum is used as the input radiation, to be reflected and 
reprocessed by an assembly of cold, optically thick clouds. 
These clouds are distributed over a spherical region from R min) 
just outside the outer boundary of the innermost region with 
the CFR cloudlets, to some maximum radius R max , which is to 
be defined shortly (see § 3). Since the lower frequency IR- 
optical portion of the primary synchro-Compton spectrum has 
been filtered out by the CFR cloudlets, the radiation pressure 
on the clouds in our region due to free-free absorption is sig- 
nificantly reduced, and consequently in this region optically 
thick clouds larger than the CFR cloudlets can exist (see 
CFR92; ST92). 

The outcoming observed spectrum I obs is calculated from the 
radiative transfer equation for an assembly of clouds iso- 
tropically distributed around the modified synchro-Compton 
source (ST92) 

r Tdd 

Ws) = U0)<T"“ + I . (6a) 

Here /,„(<)) is the input spectrum, S cld is the source function, 
and T c]d is the optical depth of the cloud system, 

t cld = ds'nn cld (s')& d {s ') , 

Jso 

where s is the path length variable and s 0 is R min . Note that T cld 
is independent of r vlcld , the optical depth of one cloud, when 
the clouds are optically thick. 

The outcoming spectrum which is due to reflection from the 
clouds is 

dR 

Rmin 

X sin 0e-^[s, dncu 

{ P /* * _ arcsin (R sin $IR mtx ) ~n Tl 

eXP L J s=o dS Sin < 6b > 

where Q is the angle subtended by the central source and the 
observer at the cloud. The parameter S ref [n th , / int (0), /] (see, 
for example, LW88) is the source function due to reflection 
from an individual cloud at R. For example, the solution of the 
radiative transfer equation for low energies is given by S ref = 
JJOXI ~ nA)/(1 + y/*)/R 2 , where e = k v /(k v + kJ. The'free- 
free opacity k v depends on the density n th m H of a cloud, which 
is derived from equations (3) and (4). The additional effects of 
the relative orientation of the cloud with respect to the central 
point source and the observer are taken into account by the 
method of phases. For example, the decrease in reflection from 
a single spherical cloud is proportional to nl 2 Jd (l + cos 6), 
similar to the solution for the luminosity of the moon phases. 
This decrease is crudely of the same order of magnitude as that 
of a circular slab tilted at an angle of 9 = 0°-60° (George & 
Fabian 1991, Fig. 12). The occultation of radiation after it was 
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reflected from the cloud is taken into account in the second 
exponential in equation (6b). 

The thermal component of the outcoming radiation due to 
absorption and reemission by the clouds is calculated from 


C Rm*x Ck ff 

dR 

dQ 

JRmln 

Jfl = 0 


H c id 2nR sin 0(i r , hermal 


x f xp [-£ 


- HTcsin (R sin 0/R m „> 


ddsin $Rn cld nl 2 ld 


}]• 


(6c) 


where thermal C^id(^)] is the thermal radiation flux from the 
surface of the clouds. 

The observed total spectrum is 


^obs fin(0)e T ' ld + / ref + / t hermai • (6d) 

This is a shorthand version of equation (6a). The integral of J obs 
gives the total observed luminosity L which should be equal to 
the integral over the primary spectrum, because we assume 
that there is no significant amount of viscosity among the 
clouds. 

The covering factor is defined as 

/(*)=[ n M nl 2 cM dR. (7) 

Then,/ = /(R max ) gives the total covering factor of the cloud 
system. The four unknowns n th , T cld , n cld , and I c)d , which deter- 
mine S ref and F therm a „ are obtained by solving equations (2), (3), 
(4), and (7). Once our free parameters L - An J dvl in ( 0), M, and 
f are specified the total spectrum is uniquely determined. For 
convenience, in our calculations an effectively equivalent 
choice of T cld (L, M,f), n cld (T cld , M), and / are adapted as our 
free parameters (see ST92 for the details). 

The integrations in equations (6b) and (6c) are solved 
numerically. The results, combined with the unabsorbed 
portion of the input spectrum, can be compared with the slab 
case where = 0.57, „ + S, d + a(7?L. - Rl in )F lh „ m „. 


3. RESULTS 

The basic equations (2), (3), (4), and (7) are solved for the 

spherically symmetric case. The results show the radial dis- 

tribution of our characteristic parameters to be T cld (R) oc 

R 1/2 > n th(R) oc R 2 , and /(R) oc log (R) to R 1 . A small cover- 

ing factor /=/(R max ) < 1 and large optical depth per cloud 

T vicid > 1 are assumed. The physical implication is that the 
luminosity is less than the Eddington luminosity. Note that for 
f (f^max) > 1 the central source is heavily attenuated and the 
radial dependence becomes more complex. The final spectrum 
was found to be almost insensitive to the radial dependence of 
n ih and/ (R) but extremely sensitive to the radial dependence of 
Jid- A positive observational check on the above results may 

e made by using CFR92. The number density of a CFR cloud- 

et near the central object is approximately n th = io 16-18 cm -3 
for an absorption column of N H < 10 21 cm' 2 for most Seyfert 

^4 ^ . W ^ e density of the broad-line region clouds at 
■ ~ V A e *l S > ~ 10 , cm 3 S « Ferland & Rees 1988). Then 

we find that n, h oc R -2 -R -2.66 ' 

Using the R dependence of the physical variables as shown 
atwye reflection is found to be most effective at intermediate 
radu. At a r^resentative radius of R = 30R, the solutions of 
equations (3), (4), and (7) give T cU ~ 3 x 10 5 * K and n th ~ 



Fig. 1. The spectra of radiation coming out of a spherical distribution of 
cool clouds are shown with the solid and dot-dashed curves, for the covering 
factor f — 0.5 and 0.9, respectively. The dotted curve shows the primary 
synchro-Compton spectrum. The chain-dashed curve is the spectrum from the 
innermost region of the central engine where very small CFR cloudlets are 
present (see the text). The dashed curve shows the spectrum coming out of a 
slab-corona system. The average temperature of the cool component is set at 
r cld = 3 x 10 s K. 


10 cm 3 * . As our free parameters the values L ~ 10 43 ergs 
s and M ~ 10 7 * * M 0 are adopted, which are representative of 
typical Seyfert nuclei such as NGC 4051 with L/2-10 
keV) = 10 41 - 5 ergs s _1 (e.g., CFR92; Ferland & Rees *1988). 
Furthermore, with /~ 0.9 the additional equation (5) gives 
hid = 10 cm and n cJd = 10 39 cm -3 . As y c i d the Keplerian 
velocity is adopted. 

As a starting point in our calculations of the outcoming 
spectra we assume that the average cloud temperature 7^ ld = 
3 x 10 5 K, the average number density of the matter in the 
cloud n th = 10 16 cm 3 , and the photon index for the synchro- 
Compton power-law spectrum T = 2. These values are chosen 
to be consistent with typical Seyfert parameters obtained from 
observation, such as L ~ 10 43 ergs s -1 and M ~ 10 7 M 0 . The 
lower and upper limits to the radial distance, R min and R max , 
are set to 10R 9 and 1000R fl , respectively. These boundaries are 
determined by the distances where the X-ray reflection ceases 
to be effective (ST92). The results are summarized in Figure 1. 
The dotted curve shows the primary synchro-Compton power- 
law spectrum. This spectrum is modified by the absorption of 
the IR-optical portion by the extremely small CFR cloudlets 
with temperature = 10 6 K. The modified spectrum is shown by 
the chain-dashed curve. Note that the absorbed IR-optical 
radiation is reprocessed and reappears as a hot thermal UV 
bump. This modified primary spectrum is used as the input to 
be reflected and reprocessed by an assembly of clouds. Since 
the clouds are optically thick in X-rays, only the unabsorbed 
portion, /| n (0)e 0 f th e original power-law source is 

observed. The spectra of the total outcoming radiation 
obtained from equation (6d) is shown by the solid curve for the 
total covering factor/ = /(R m J = 0.5, and by the chain-dotted 
curve for/= 0.9. For comparison the dashed curve shows the 
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Fig. 2. — Here the 5-10 keV section of the spectrum in Fig. 1 is enlarged. 
The upper dotted curve is now for/- 0.01. Otherwise the notation is the same 
as in Fig. 1. 


spectrum of radiation reflected from a single cold slab, using 
the same input spectrum (LW88). Figure 1 is enlarged in 
Figure 2 to show the details of the 5-10 keV X-ray region. The 
upper dotted curve is now the spectrum for / = 0.01. Otherwise 
the notation is the same as in Figure 1. The photoelectric effect 
is represented by the Fe 6.4 keV emission line and 7.1 keV 
edge, for demonstrative purposes. 

In Figure 1 the broadening of the UV bump extending from 
the UV to the soft X-ray excess region, as well as the effect of 
free-free absorption and Compton recoil, are clearly seen. In 
Figure 2 the “flattening” of the power-law spectrum from 
T = 2 to F = 1.7 + ®;7 in the X-ray band (e.g., Zdziarski et al. 
1991) is evident. Also noted is the increase in the equivalent 
width with increasing/. Note that the equivalent width for the 
cloud model with f~ 0.7 is similar to that of the slab case with 
/~ 0.5, while the Fe edge is much smaller for the cloud model. 
Even without the washing-out effect of the Doppler shifts, etc., 
the edge is not noticeable as/is decreased below ~0.5. 

Previous authors have found good agreement with the X-ray 
spectral data of some AGNs when the reflection from a single 
slab is adopted (e.g., LW88; George & Fabian 1991, and refer- 
ences therein). Qualitatively, our model also gives similar 
X-ray spectra when 0.1 </< 1. However, our cloud model is 
more flexible than the disk model because it generally covers a 
wider range of the covering factor. For instance, for the stan- 
dard disk model /is fixed at 0.5. Quantitatively, with the same 
value of/ the major spectral effects of the cloud configuration, 
as compared with the slab geometry, are a significant increase 
of the UV bump and a significant decrease of the reflection 
features, such as the Fe line-edge system and the hard X-ray 
bump. The latter effect is caused by the obscuration of the 
reflected component by the clouds between the reflecting cloud 
and the observer. 

Most of the covered portion of the input spectrum is 
absorbed by the clouds and reradiated as thermal radiation 
producing a large UV bump. Therefore, we calculated the total 
integrated bolometric luminosity for the most extensively 


studied rapidly variable Seyferts, NGC 4051 and NGC 6814, 
in order to check whether it would exceed the Eddington lumi- 
nosity. It is found to be less than or roughly equal to the 
Eddington limit. In this calculation the X-ray portion of the 
spectrum was set equal to the observed X-ray intensity of these 
objects. Also assumed are a large /, e.g., ~0.9, and the black 
hole mass of 10 6 — 10 7 M 0 . 

4. DISCUSSION AND CONCLUDING REMARKS 

Let us compare the variability behavior of the cloud model 
and the thin disk models. Intensity variations of the broad 
thermal UV bump produced by the clouds are temperature- 
and frequency-dependent. For the cloud model a “ turnoff” of 
the central X-ray emitting source will result in cooling of the 
clouds. The time scale is derived from dE/dt = -Anl^ d a Q T* ]d , 
where E is the total energy of one cloud and T cld (E) = 
2>E/{Anll Xd n th k B ) is found from the ideal gas law. Solving this 
equation we get 

3 x 10 16 

t c >— 3 s, (9) 

1 cld 

for Seyfert nuclei, e.g., NGC 4051 (ST92). For instance, for very 
soft X-rays at 7^ ld = 10 6 K t c ~ 0.03 s, while it is nearly a year 
for the UV radiation at T cId = 1000 K. Therefore, with our 
model, while the high-energy end of the bump should vary 
almost simultaneously with the harder X-rays, similar rapid 
variations are not expected at the low-energy end of the bump 
lying in the UV band. We carefully examined the spectral 
behavior of the NGC 4051 data near the B band (Done et al. 

1990) which shows no rapid variability during the X-ray flares, 
and find that the low-energy end of the UV bump already 
appears there. Therefore the lack of rapid variability at this 
wavelength is consistent with our cloud model. In our model 
the power-law portion of the IR-optical radiation is due to 
different mechanisms (Tsuruta 1990), and hence no correlation 
is predicted between the IR-optical-UV bands and the X-ray 
band. On the other hand, definite correlated variations are 
predicted among the IR-optical-UV and X-ray bands when a 
variety of disk models are adopted; for instance, the standard 
thin-disk model (Shakura & Sunyaev 1977) and a disk with a 
hot spot (Abramowicz 1990). Simultaneous X-ray continuum 
and Fe line variations (within ~300 s) observed in NGC 4051 
(Kunieda et al. 1990) are consistent with our cloud model if the 
line is produced at <30R 9 , because then, with M ~ 10 6 M 0 
predicted for NGC 6814, the delay in the line variation should 
be < 300 s. 

The cloud model might also be applied to galactic X-ray 
sources which are black hole candidates, because the under- 
lying physical mechanisms governing these objects are con- 
sidered to be similar. In the case of the galactic black hole 
candidate Cygnus X-l, our cloud model agrees better with 
observations than disk models, because the weaker and sym- 
metrical Fe line feature predicted from our models with/< 0.2 
is consistent with the observed data of this source, while the 
slab reflection models predict this feature to be stronger than 
and different in shape from the observed Fe line (Done et al. 

1991) . 

Another means to distinguish between the disk and cloud 
models comes from the observation of polarization. In thin 
disks scattering atmospheres should induce a high degree of 
polarization (Netzer 1990, and references therein). On the other 
hand, our cloud model predicts very low polarization, because 
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